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We discuss in detail the derivation of the leading four-nucleon force in chiral effective field theory 
using the method of unitary transformation. The resulting four-nucleon force is given in both 
momentum and configuration space. It does not contain any unknown parameters and can be used 
- . . in few- and many-nucleon studies. 
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I. INTRODUCTION 



Chiral effective field theory (EFT) is a powerful tool for analyzing the properties of hadronic systems at low energies in 
a systematic and model-independent way. It exploits the approximate and spontaneously broken chiral symmetry of 
QCD which governs low-energy hadron structure and dynamics. Over the past years, considerable progress has been 
achieved in understanding the structure of the nuclear force in this framework, see [l|, 01 for recent review articles. In 
particular, two-nucleon force (2NF) has been worked out and applied in few-nucleon studies upto next-to-next-to-next- 
to-leading (N'^LO) in the chiral expansion [1, Three- nucleon force is currently available at next-to-next-to- leading 
order (N^LO) [1, see 0, H, @ for some examples of recent few-nucleon studies at N^LO. Although the results for 
most of the investigated few-nucleon observables look promising, it is necessary to go to N'^LO in order to test the 
convergence of the chiral expansion. Increasing the chiral order is also expected to bring new insights into certain 
■ existing puzzles in the 3N continuum such as e.g. the large deviations from the data for the differential cross section 
ly-j in some deuteron breakup configurations, see [lO| for a recent work on this subject. The full N"^LO analysis requires 
^sj . the incorporation of the subleading contributions to the 3NF which are currently being worked out. In addition, 
' one has to take into account four-nucleon force (4NF) which starts to contribute at this order. In Ref. [ll| we 
, already presented the expressions for the leading 4NF. It is governed by the exchange of pions and the lowest-order 
y—^ ' nucleon-nucleon contact interaction and includes effects due to the nonlinear pion-nucleon couplings and the pion 
t — self-interactions constrained by the chiral symmetry of QCD. The obtained 4NF is local and does not contain any 
unknown parameters. In this work, we describe in detail the derivation of the 4NF using the method of unitary 
^ transformation (T^ . ITsll. fo r which a new formulation is presented. This new formulation is considerably simpler than 
the one given in Refs. [l2l[l3j and much more convenient for practical applications. We furthermore study the effects 
of the additional unitary transformations which affect the form of the effective Hamilton operator at N'^LO. The 
appearance of such transformations is a new feature at this order in the chiral expansion. 



Our manuscript is organized as follows. In sectio n HIl we discuss the chiral power counting in the few-nucleon sector 
following the original formulation by Weinberg 

El 

[T5| . We express the power of the low-momentum scale in a 
different form as compared to Refs. [3, fisj which is particularly useful for applications based on algebraic rather 
than diagrammatic methods. Combining this new formulation with the method of unitary transformation we obtain 
the formal expression for the effective nuclear Hamilton operator in a compact form. In section Hill we use this new 
formulation to derive the N'^LO contributions to the effective nuclear Hamiltonian and evaluate the corresponding 
four-nucleon operators which give rise to the 4NF. We employ a large class of additional unitary transformations 
acting on the nucleon subspace of the full Fock space and investigate their impact on 3NFs and 4NFs. We end with 
the summary and outlook. Finally, appendix [X] contains some details related to the N^LO contributions involving 
subleading vertices while the expressions for the 4NF in configuration space are given in appendix [Bl 
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II. NUCLEAR FORCES IN CHIRAL EFT USING THE METHOD OF UNITARY TRANSFORMATION 



A. Chiral power counting 



Power counting is a crucial ingredient of any effective field theory which allows to organize various contributions to the 
scattering amplitude according to their relevance. Throughout this work, we adopt Weinberg's power counting which 
is based on naive dimensional analysis^. Traditionally, chiral power counting is formulated in terms of topological 
quantities such as the number of loops and, in the few-nucleon sector, the number of participating nucleons and 
disconnected pieces [l6l |. This allows to immediately identify all relevant diagrams at a given order and 

appears to be particularly useful for calculations based on the Feynman graph technique. For algebraic approaches 
such as the method of unitary transformation (l^ . [l3| it is more convenient to express the power counting in a different 
form which will be given below. For the sake of completeness, let us first recall the derivation of Weinberg's power 
counting result following his original work [T^ . ITsj . Consider an arbitrary A^-nucleon irreducible time-ordered diagram 
(i.e. the one that does not contain iterative contributions) without external pions. For the sake of generality, we 
allow for cases in which N out of N nucleons do not interact. To obtain the corresponding power ly of the soft scale 
Q one has to count the number of the independent three-dimensional momentum-space integrations / — Vi with 
/ (Vi) being the number of internal lines (vertices of type i), the number di of derivatives or pion mass insertions at 
each vertex, the number D of energy denominators and the number of phase-space factors associated with exchanged 
pions: 

ly = 31 - 3 (^V, - 1^ - D + {d^ - j) -3iV. (2.1) 

Here, pi denotes the number of pion fields at a vertex of type i. Notice that we follow the standard convention to 
define the potential. In particular, we do not count the overall (5-function that enters the definition of the S-matrix 
element but do count the additional (5-functions associated with the noninteracting nucleons. To bring Eq. (|2.ip into 
the standard form one can use the following topological identities. First, the number of intermediate states D can be 
expressed as: 



Secondly, the number of loops L is given by 



= ^ - 1 . (2.2) 



L = I-Y,y^ + G-N, (2.3) 



where C refers to the number of separately connected pieces. Notice that each noninteracting nucleon is regarded as 
a separately connected piece. The last identity we need reads 

2/ + 2iV = ^y,(p, +nO + 2iV. (2.4) 

i 

Here, Ui is the number of nucleon field operators at a vertex of type i. Using Eqs. p.2p . (|2.3|) and (|2.4|) . one can 



■ using _tjcis. 

rewrite Eq. (|2.ip in the form that appears in Weinberg's original papers [iSl. [T7|: 

i/ = 4-7V + 2(L-C)-f ^1^,A,, A, = + -2. (2.5) 

i 

There is one subtlety in the above expression which needs to be addressed: according to Eq. (|2.5p . the chiral dimension 
v for a given process depends on the total number of nucleons in the system. For example, one-pion exchange in the 
two-nucleon system corresponds to iV = 2, i = 0, C = 1 and Aj = and, therefore, contributes at order v = 0. 
On the other hand, the same process in the presence of the third (spectator) nucleon leads, according to Eq. ()2.5p . to 



Notice that alternative counting schemes for contact two-nucleon interactions are currently being explored. 
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!/ = — 3 since iV = 3 and C = 2. The origin of this discrepancy hes obviously in the different normahzation of the 2N 
and 3N states: 

2A^ : {pi P2\piP2) = 5^{pi - Pi ) S^{p2' - P2 ) , 

3iV: (piP2P3|piV2'p3') =<5'(Pi'-Pi)<5'(p2'-P2)<53(p3'-P3). (2.6) 

It can be circumvented by assigning a chiral dimension to the transition operator rather than to its matrix elements 
in the A^-nucleon space. Adding the factor 3iV — 6 to the right-hand side of Eq. (|2.5p in order to account for the 
normalization of the 7V-nucleon states and to ensure that the LO contribution to the nuclear force appears at order 
1/ = we obtain 

j/ = -2 + 27V + 2(L-C)+^F,A, . (2.7) 

i 

We will now rewrite this expression in a different form which is better suited for the method of unitary transformation. 
To that aim we combine Eqs. (|2.3p and (|2.4p and express 2(L — C) as 

2(L -C) = -2N + ^ V. ip^ + n, - 2) . (2.8) 

Substituting this expression into Eq. (|2.7p leads to 

j/ = -2 + ^ViK, , = rfi + -n, +pi - 4. (2.9) 

Clearly, the quantity Ki which enters this expression is just the canonical field dimension of a vertex of type i (up to 
the additional constant —4) and gives the inverse mass dimension of the corresponding coupling constant. It should 
further be emphasized that Eq. (|2.9p can be obtained immediately by counting inverse powers of the hard scale A 
rather than powers of the soft scale Q (which is, of course, completely equivalent). Indeed, since the only way for the 
hard scale to be generated is through the values of the low-energy constants (LECs), the power v is just the negative 
of the overall mass dimension of all LECs. The additional factor —2 in Eq. (|2.9p is a convention to ensure that 
the contributions to the nuclear force start aX v = Q. While Eq. (|2.9p does not say much about the topology and is, 
therefore, not particularly useful to deal with diagrams, it is very convenient for algebraical calculations. In particular, 
it allows to (formally) reduce the chiral expansion to the ordinary expansion in powers of the coupling constant. The 
role of the coupling constant is played by the ratio Q /A, and the power of the coupling constant for a vertex of type i 
is given by Ki. Clearly, for perturbation theory to be applicable it is necessary, that only nonrenormalizable (i.e. the 
ones with Ki > 0) interactions appear in the Lagrangian. This is guaranteed by the spontaneously broken chiral 
symmetry. 



B. Application to the method of unitary transformation 

We are now in the position to apply Eq. (|2.9p to derive nuclear forces from the effective chiral Lagrangian using the 
method of unitary transformation. The starting point is the time-independent Schrodinger equation for interacting 
pions and nucleons 

{H^ + Hi)\^)=E\^), (2.10) 

where Y^) denotes an eigenstate of the Hamiltonian H with the eigenvalue E. Let ?/ (A) be projection operators onto 
the purely nucleonic (the remaining) part of the Fock space satisfying rf' = rj, = X, rjX — Xi] ~ and A + 77 = 1. 
To study nuclear systems below the pion production threshold it is advantageous to project Eq. (|2.10p onto the ?/- 
subspace of the full Fock space. The resulting effective equation can then be solved using the standard methods of 
few- or many-body physics. This reduction can be achieved via an appropriately chosen unitary transformation 

Following Okubo fl^, the unitary operator U can be parametrized as 



_ / 77(l + ^t^)-i/2 -At(l + AAt)-i/2 \ 
~ A(l + At^)-i/2 A(l + AAt)-i/2 J 



(2.12) 
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with the operator A = XArj. The operator A has to satisfy the decouphng equation 

X{H- [A, H] - AHA) r]^0 



(2.13) 



in order for the transformed Hamiltonian H to be of a block-diagonal form. The effective ry-space potential V can be 
expressed in terms of the operator A as: 



V^7^{H- Ho) = V 



(1 + + A^H + HA + A^HA){1 + A^A)-^/'^ - Hq 



1] 



(2.14) 



The potential V can be derived perturbatively based on the chiral power counting. In Refs. [13, ll3| this was achieved 
using Eq. (|2.7p via a two-dimensional recursive process which required rather tedious calculations in the intermediate 
steps. As an alternative, one can use Eq. (|2.9p to express Hj as 



Hi = Y,H'' 

K = l 

with K being defined in Eq. (|2.9p . Assuming the following expansion for the operator A 

oo 

one immediately obtains from Eq. (|2.13p : 

I K — 1 1 K — 2 K~j~l 

A'' = - —xl^H'' + Y,H'A''-'-^A''-'W-Y^ A'WA""- 



(2.15) 



(2.16) 



(2.17) 



1=1 j=i 



Here, Erf {E\) refers to the free energy of nucleons (nucleons and pions) in the state ry (A). The expression for the 
effective potential follows immediately by substituting Eqs. (|2.15p and (|2.17p into Eq. (|2.14p . 



III. DERIVATION OF THE LEADING FOUR-NUCLEON FORCE 



A. Effective Lagrangian 

The effective chiral Lagrangian for pions and nucleons has the form, see e.g. [19|, |2 

C^^ = ^ [d^Ud^U^ + Ml{U + C/t)] + (3.18) 
C^N = (^Do -^a-u)N+... , 

Cnn = -'^Cs{N^N){N^N) ~ ]^CT{N^aN) ■ {N'^aN) + . . . , 

where only those terms are shown explicitly which contribute to the leading 4NF. Here, = 92.4 MeV [gA = 1.267) 
is the pion decay constant (the nucleon axial- vector coupling), N represents a non-relativistic nucleon field and a 
denote the Pauli spin matrices. The low-energy constants Cs and Ct determine the strength of the leading NN 



short-range interaction [ij, [l^- Further, the SU(2) matrix U — v? collects the pion fields, D^' = + \ [u\ 9'" 



u\ 



denotes the covariant derivative of the nucleon field and — iu^d^JJu^ . The first terms in the expansion of the 
matrix [/(tt) in powers of the pion fields take the form 

Uliz) = 1 + . - AtTt^ - ^(t • 7r)3 + ^^^71-4 + . . . , (3.19) 



where r denote the Pauli isospin matrices and a is an arbitrary constant. Notice that only the coefficients in front 
of the linear and quadratic terms in the pion field are fixed uniquely from the unitarity condition C/^[/ = 1 and the 
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proper normalization of the pion kinetic energy. The exphcit a-dependence of the matrix U represents the freedom 
in the definition of the pion field. Clearly, all measurable quantities are be a- independent. 

Expanding the terms in Eq. (|3.18p in powers of pion fields and applying the canonical formalism leads to the following 
interaction terms in the Hamilton density: 



= 

= 



4a- 1 



F2 



8a -1 



+ i Cs{N^N){N^N) + ^Ct [N^aN] ■ (N^aN) , 



1 



2a 



(t • tt){tz(t ■ Vtt) + a7r^(x(T • Vtt) 



N. 



TVV • (tt X TV)N 



[N\t X 7r)iV] • [N\t X 7r)iV] . 



(3.20) 



Here, the superscripts of Ti. refer to the canonical field dimension k defined in Eq. (|2.9p . Again, we only show those 
terms in the Hamilton density which are relevant for the present calculation. Notice that the last term in Eq. p.20p 
is absent in the effective Lagrangian in Eq. p.lSp and arises through the application of the canonical formalism, see 
e.g. [13. 



B. Orders v < A. 



The structure of the effective Hamilton operator at leading order can be worked out straightforwardly applying the 
projection formalism as described in section [HB] to the lowest-order (i.e. with = 0) Hamilton density in Eq. (|3.20p . 
This leads to 



ff2 _ rrl ttI 



(3.21) 



Here, the superscript of A refers to the number of pions in the corresponding intermediate state, E-^ is the total energy 
of n pions in the state A", — X^iLi '^i where Ui = \/ qi'^ + is the energy of the pion with momentum g^. The 
subscripts a and b of H^f^ denote the number of nucleon and pion field operators, respectively. They are introduced 
in order to clarify the meaning of various terms in the effective potential. Clearly, the terms in Eq. (j3.2ip only give 
rise to one- and two-nucleon operators. The first corrections to the effective Hamiltonian arise at order v = 2. The 
relevant terms at this order read: 



Xl \2 \1 

ttI ffl Wl 



Ett Et^ 



Ai 

E-r; 



^21 -f^H^^Q 



Ai_ 

Ett 



E. 
1 



21 



1 A^ 

2^21 ~^^21 V -^^21 



1 A^ 



(3.22) 



Here we onl y li st those terms which can generate four-nucleon operators. The complete list of terms at this order can 
be found in [l2l . [2]| . The contribution to the 4NF from the operators in Eq. (|3.22p can be represented schematically by 
disconnected diagrams shown in Fig.[TJ It should be understood that diagrams in the method of unitary transformation 
have a different meaning compared to the ones arising in the context of time-ordered perturbation theory and serve 
merely to visualize the topology corresponding to a given sequence of operators H'^^^. We further emphasize that 
diagrams shown in this work represent the sum over all possible "time orderings" as depicted in Fig. O 

The contributions to the 4NF from diagrams in Fig.[T]have been considered in Ref. [5] using time-ordered perturbation 
theory. In this approach, only the first and the fourth terms in Eq. p.22p contribute to the effective Hamiltonian. 
Evaluating matrix elements of these terms corresponding to diagrams shown in Fig. [2] one obtains non-vanishing 
contributions to the 4NF. As pointed out in Ref. Q, the resulting 4NF cancels against the recoil corrections to the 
2NF upon the iteration in the dynamical equation. The same sort of cancellation also takes place for the 3NF at 

v=2 sm. 
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FIG. 1: Disconnected 4N diagrams at order v = 
dots denote vertices with Ai = from Eq. (|3.20p . 



2. Solid and dashed hnes represent nucleons and pions, respectively. Solid 



+ 



+ 



FIG. 2: Different "time orderings' 
nucleon lines are not shown. 



for the first diagram in Fig. [T] respectively. Graphs resulting by the interchange of the 



In the method of unitary transformation, one has to take into account all terms in Eq. (|3.22p . It is easy to verify 
that the resulting 4NF vanishes due to an exact cancellation between the different terms in this equation. A similar 
cancellation occurs also for the 3NF [13, H^l, see also [l^l, and for disconnected 2NF diagrams [2l| at this chiral order. 
The corrections to the effective Hamilton operator at order = 3 within the method of unitary transformation are 
discussed in detail in [l3l | . They give rise to at most tree-nucleon operators and will not be discussed in this work. 



C. Order = 4. 



The first non-vanishing contributions to the 4NF arise at order i/ = 4. For the sake of a better overview, it is useful 
to divide the total contribution at this order into pieces with the same combinations of the coupling constants and 
Cs,T- This leads to eight classes of terms which are discussed in detail below. 



Class-I contributions proportional to g\. 

This class of contributions arises from all possible 4NF diagrams involving six vertices from the first line of 
Eq. ((3:201) . Using Eqs. ((2l4)) - ((2l7)) one obtains: 



A^ A^ A^ A^ A^ 

H2iTrH2ir]H2i—H2i—H2i—2H2i + H2i—H2i'qH2i—H2i—2H2i—H2i 

y\ y\ y2 y\ y\ y\ y2 y\ 

+ H2i-^H2iriH2i—2H2i—H2i — H2i + ^^2i^-f^2i ^21-^-^^21 ^21^ -^^21 

J^n J^TT -C/TT J^TT -C'tt 

^1 I 
- H21 H^i vHli— H^i 1] H^i - ^ -^H^iV H21 H^i V H^i 

~ 4 -^21 -^^21^ H21 H21 T] H21 H21 - Hl^ — — — H21 — H21 — H21 



Tjl '^^ Tjl TTl TTl Tjl Tjl 

-"21 p ^21 p -"21 p -"21 p ^21 p -"21 

J^TT -L^TT -C^TT -L^TT -L^TT 



77 + h.c. 



(3.23) 



The corresponding contributions to the 4NF can, in principle, be evaluated straightforwardly by calculating 
matrix elements of the operators in Eq. (|3.23|) for all possible "time orderings" of diagrams shown in Fig. [3l 
Before giving the explicit results, it is important to address the issue whether the contribution V^^^ is defined 



7 



FIG. 3: Class-I contributions to the 4NF. For notation, see Figs. [T] and [2] 



unambiguously. Generally, an effective Hamilton operator is defined modulo unitary transformations.^ It 
should, therefore, be understood that terms in Eq. (j3.23p correspond to one particular choice of the unitary 
transformation U . Indeed, Eq. (|2.12p does not yield the most general parametrization of the operator U . The 
resulting effective Hamilton operator rjHrj can be further modified via subsequent unitary transformations acting 
on the 77-space. To be specific, consider the unitary transformation of the form: 



?7 = e-^ , 

where S is an anti-hermitian operator acting on the 77-space, S = rjSrj, 
S — aiSi + 0282 with ai^2 being arbitrary real numbers and 



(3.24) 

-S. Further, let S be given by 



Si = T] 

82=1] 



H2i^H2ir]H2i-^H2i - H2i^H2iri H21—H21 77, 

I A^ A^ A^ A^ 

H2i-^H2i—H2i^H2i - H2i-^H2i—H2i—H2i 



(3.25) 



The operators Si and ^2 are the only possible ones that are invariant under the time-reversal operation and 
can be constructed out of four vertices i?2i provided the iJ^r's are only allowed to appear in the denominators. 
Acting with the transformation U onto the lowest-order effective Hamilton operator. 



1 1 



(3.26) 



where fffp is the nonrelativistic kinetic energy term, the following additional terms in V^^'^ are generated: 



-ai rj 

- a2r] 
f h.c. 



A-^ A^ A-"^ A-^ A-^ A-^ 

H21 H21 rj H21 -^Hl-^V H21 -^^21 - H21 — H21 77 i?] 1 H21 rj H21 — H21 

\i A^ A^ A^ A^ A^ A^ 

Hj-rr rjTT rjTT tLt tliTT hit. JliTT rjTT 



(3.27) 



Here, the ellipses refer to terms involving H2Q (fffo) which give rise to class- VIII (class-IV) corrections and 
will be considered below. The class-I contributions to the effective Hamiltonian, therefore, seem to be defined 
modulo the arbitrary constants ai and a2- We will now demonstrate that there exist one particular choice for 



^ Other sources of ambiguities related to field redefinitions in the Lagrangian and to the choice of the dynamical equation are discussed 
for one- and two-pion exchange potentials in Ref. |25|| . 
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(a) 



(b) 



(c) 



(d) 



(e) 



(f) 



(g) 



FIG. 4: Selected three-pion exchange contributions to the 3NF at order u — 4 (graphs (a)-(f)) and one-pion exchange 3NF at 
order i/ = 4. Filled circle denotes the leading ttNNNN vertex with Ai = 1. For remaining notation, see Figs. [l]and[2l 



these constants which is strongly preferable. To that aim consider three-nucleon forces generated by terms in 
Eqs. p.23|) . (|3.27p and involving one-pion exchange between two nucleons, see Fig. |4l diagrams (a)-(f). We can 
generally express the structure of the 3NF in the form 



^3N — ^'^3N ^3N ■ 



(3.28) 



Here, Mgjy represents the spin, isospin and momentum structure which results from vertices entering a diagram i 
and ^3*^ denotes the corresponding energy denominators. Evaluating matrix elements of operators in Eqs. (|3.23p 
and l|3.27p for all possible "time orderings" of diagrams (a)-(f) in Fig. 0] yields 



^3N 
^3N 



3N 



3N 



3N 



4 



^3N — 



3N 



{a,ij}^{c,J,i} 

- 2a2 1 + 2ai 

2 2" 



1 + 2q;i 



1 

■f 

1 - 2^2 



\-2ai ^l-2ai 



8a2 



1 + 2ai 



1 

-f — 



1 + 2q!i 1 + 2ai 



l + 2ai 



1 + 2ai 



, 1 + 2ai 

OJiUJ^ (Ujj +UJk) " i^f^j i^^i + t^k ) 



l + 2ai+4a2 „1 



3 

2ai + 4a2 



3N 



{d,i,j}^{f,j,i} 



(3.29) 



The contributions from crossed-box diagrams (a-c) exhibit the expected dependence on the pion energy ujk which 
leads to the usual static one-pion exchange potential (OPEP) c>c — l/{qk^ + M^). The additional factor of 
uj^^ results from the pion phase-space factors. Notice further that V^'^''^ do not depend on the parameters ai,2'^- 
The contributions from the box diagrams (d-f) do, however, depend on ai.2 and involve terms that cannot be 
associated with the static OPEP between the second and the third nucleons. For example, terms proportional 
to WjT^ in the last three lines of Eq. (|3.29p lead to the 3N potential proportional to (qk ^ + M^)^^^^ which has 
a cut starting at = —M^. The presence of such a cut which is not produced by multi-particle intermediate 
states appears, at first sight, rather puzzling. One should, however, keep in mind that we are dealing here with 
the non-iterative part of the scattering amplitude which gives rise to nuclear forces. Its analytic structure needs, 
in principle, not to coincide with the one of the full amplitude. A more serious problem due to the appearance of 
the "unphysical" terms in Eq. (|3.29|) is non-renormalizability of the corresponding 3NF. Indeed, one can easily 
verify using dimensional arguments that the loop integrals entering, for example, the oj^T'^ i'^k'^) t^rms give rise 



^ The independence on cJi follows directly from Eq. I|3.27| l. 
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to cubic and linear (quadratic and logarithmic) ultraviolet divergences. The only relevant counter term which 
is available at this order and which can be used to absorb the infinities is H^, with one pion and four nucleon 
field operators. It is proportional to the LEG D and enters the 3NF at order v = 3, see diagram (g) in Fig. lU 
which arises from 

- ~V —^21 + ^21 — ^4il V , (3.30) 

and leads to V^j^ — —2/uji~ The contributions to the 3NF in Eq. (|3.29p with a different dependence on 

Uk can, obviously, not be renormalized in this way. Again, it should be understood that the above mentioned 
difficulty does not affect S-matrix elements. All "problematic" divergences that enter the 3NF have to cancel 
against the divergences arising from iterative contributions to the amplitude in such a way that the resulting 
S-matrix is renormalizable in the usual sense. For the purpose of describing the few-nucleon dynamics based 
on the Schrodinger equation, it is, however, desirable to have a well-defined effective Hamilton operator. The 
difhculty with the non-renormalizability of the 3NF in Eq. (|3.29[) can be avoided if one requires that only the 
^ terms are present, i.e. the OPEP factorizcs out in diagrams (d-f) in Fig. [H This can be achieved via a 



suitable choice of the r/-space unitary transformation by setting 



2 

This choice leads to the following remarkably simple expressions: 



ai = -a2 = -\ ■ (3.31) 



3^ ~ c.^f^fe 



AT — 



— ~ ! 3 



^3N — ^3N 



(3.32) 



These considerations may remind one of the recent findings in the context of large- A^c QCD [2^ |27|, |28| . 
There it was found that the multiple-meson exchange potential derived in the energy-dependent formulation is 
inconsistent with large- A^c counting rules. The consistency could be maintained using a different (but completely 
equivalent) form of the potential based on the energy-independent formalism, see Ref. [2^ for more details. 

Choosing the parameters ai_2 as described above, it is a straightforward exercise to calculate the contributions 
to the 4NF associated with the diagrams in Fig. [31 We found that all disconnected graphs lead to vanishing 
contributions regardless of the values of ai,2- From the connected diagrams in the second row in Fig. [Sj only 
the first two generate non-vanishing contributions to the 4NF which take the form: 

I'Class-I - 



(2F.)6 + ^^2] [g.^2 + ^^2]2 [^2 + ^^2] 

(ti ■ T4T2 ■ T3 - Ti ■ X3 T2 • T4) qi ■ qi2 ■ qi2 + Ti X T2 ■ T4 qi ■ qx2 qi2 x qA ■ ^3 



Tl X T3 • T4 qi ■ qi2 qi x qi2 • (72 + Ti • 1-4 qi2 x qi ■ <T2 qi2 x 94 • 0-3 



all perm.. (3.33) 



Here, the subscripts refer to the nucleon labels and qi — pi' — Pi with pi ' and pi being the final and initial 
momenta of the nucleon i. Further, qi2 = qi + q2 = —Q3 — Qi = —Qsa is the momentum transfer between the 
nucleon pairs 12 and 34. We have verified that the obtained expression for the class-I 4NF remains unchanged 
if one considers a larger class of unitary transformations in Eq. (|3.24p with the operator S involving terms 
with pion energies in the numerators. Although the renormalizability condition does not completely fix the 
corresponding 77-space unitary transformations in that case, the remaining ambiguity does not affect nuclear 
forces at the considered order. Notice further that the contribution from the last graph in Fig. [3] only vanishes 
if the constants ai^2 are chosen according to Eq. (j3.3ip . Finally, it should be emphasized that the ambiguity 
due to the additional ?7-space unitary transformations does not show up at lower orders in the chiral expansion 
since it is not possible to construct non- vanishing anti-hermitian operators S with just two vertices H^i- 

• Class-II contributions proportional to g\. 

These contributions arise from 4N diagrams involving four vertices H21 and one insertion of the Weinberg- 
Tomozawa vertex H22, see Fig. [SJ In addition, there are diagrams with one insertion of the nonlinear pion 
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FIG. 5: Class-II contributions to the 4NF. For notation, see Figs. [T] and (2] 



and pion-nucleon interactions i?o4 and H23, respectively. Using Eqs. (|2.14p - (|2.17p one obtains the following 
contributions to the class-II effective Hamilton operator: 



1/(4) 



V 



1 A^ A^ 

2 ^2 1 Hli V 2 1 ^22 



A 



Ai 



A2 



22 K 21^2-21 



1.1 



A 



'22 



El 



■E. 



Ett 



1.1 A^ 



n A^ ^1 A 



7;Hli -f^Hl^ V Hli -ftHIi—HI^ - -HI^—^HI^ 77 HI^—HI^—HI 



1.1 A^ 



■El' 
n A^ 



1 -y 2 A^ 

2^21^^217] H22^H2i—H2i 



2—21 ^2 21 '/ -'-'21 ^ 22 -'-'21 

1 1 A"'^ 1 A"^ n A'^ 1 A"'^ T 

0-"21 77 -"21 77 -"22 77 21 77 -"21 

Z ii^r -C/TT ^TT J^TT 



A 



A2 



+ -^21 -7^E[li -t^HI 



Ai 



-'-^TT J-^TT J-^TV ^TT 



-H 



A2 



21" 



-H 



22 



H 



A 



21 



-Hn 



J-^TT ^TT ^TT ^TT 



Ai 



21" 



A 



22" 



-Hi 



21 



frl A^ rrl A^ l A^ , A^ 2 



Et, 



' Ett 



ttI A^ tti ^ A^ 2 A^ 1 



Xl X2 \3 



ttI rrl tjl tr3 rrl rrl rr3 rrl 

J 1 _ — i i r) 1 _ — J 1 Q 1 _ — J / Q 7 7 1 _ — J 7 Q 1 _ — J 7 Q rj ~ ' — 7 7 q ■ 



1 IJl 



■E^ 



Ett 



' Ett 



Ett 



Ett 



' Ett 



' Err 



+ Hl 



21 



-H 



^TT ^IX ^TT ^TT 



21" 



X4 

tfl W2 

-n^i—n, 



04 



r/l A^ rrl A^ ^ A^ 3 A^ 1 
21 " J2 01 „ ri „ n/i „ ^ ^ ') 



1 fjl A rrl A 2 A 1 A , 

0-"21 T^-"21 -"o4i=r-"2i i=r-"2i 

Z il/7r £/7r ^TT 



•q + h.c. 



'21 



^04 



E^ 



21 



(3.34) 



Similar to the previously considered class-I contributions, we employ additional 77-space unitary transformations 
with the operator S in Eq. p.24p given by 5 = ass's -f 014^8/^ -f a^^S^ where are real constants and 



S3 


= V 


Hli 




Ai 

Eti 




- -f^2\4r- 


-^^22 


Ai 




V, 


Si 


= V 


HI2 


)2 


Ai 
E, 




- -f^2iTr 


-f^2\ 


A2 


HI2 


V, 


S5 


= V 


H22 


\2 

-prHli 


Ai 


hIi 




-f^2\ 


A2 

Et! 


H22 


V- 



(3.35) 

The operators S3, S4 and S5 are the only time-reversal invariant anti-hermitian operators that can be constructed 
out of two vertices iJ2i and one Weinberg- Tomozawa vertex H22 with E-r^'s appearing only in the denominators. 
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(a) 



(b) 



(c) 



(d) 



FIG. 6: Class-II three-nucleon diagrams which provide constraints on 03, 04 and as as explained in the text. For notation, 
see Figs. [2 and [21 



The corresponding Ofi-dependent class-II contributions to the effective Hamilton operator read: 



^21 — -^21 ^ -^^21 ^-^22 



a4 T] 



as?? 



2 Tjl 

Ai 



A^ A^ A^ 

■^^21 — -f^21 V H21 —H22 



El 



21 



Ett 



H21 



Hi 



Ett 

Ai 



" Ett 

Ai 

" Ett 



A2 



A^ A^ A^ 

H 21^ H 21 Tj H 22-^ H 21 H21 



21 ^2lV ^21 ^21 j^2 ^22 

A^ A^ A^ 

^21^ H21 rj H2i^H2i — H22 



(3.36) 



77 



h.c. 



where the ellipses refer to terms involving an insertion of either the nucleon kinetic energy i?|g or the contact 
interaction H^q (class-V contributions). The constants as, 0:4 and as are constrained by the requirement that 
the OPEP factorizes out in the 3NF diagrams shown in Fig. [51 This guarantees that all ultraviolet divergences 
arising in the corresponding loop integrals are absorbable into a redefinition of the LEG D entering the 3NF (g) 
in Fig.m We emphasize that the 3NF diagrams shown in Fig. [S] are not the only possible class-II diagrams with 
a single pion being exchanged between the first two and the third nucleons. All other 3NF contributions which 
are not shown in Fig. [6] are found to be proportional to w^T^ irrespectively on the values of ai and, therefore, 
do not constrain these constants. Evaluating matrix elements of the operators in Eqs. (|3.34p and (|3.36p for all 
possible "time orderings" of diagrams in Fig. [H] yields: 



V-. 



3N 



LUiLUjLUk 



3^ 



{a,iJ}^{b,jA} 



V-. 



I "3 + as 1 — 2a4 



3N 



4a5 2 3 — 2a4 + 4a5 



, ^3 + as 



U!iUJjUJk 



3Ar 



3N 



{c,i,j}^{d,j,i} 



(3.37) 



In order to simplify the above expressions, we have taken into account the pion kinetic energy arising from the 
time derivative in the Weinberg- Tomozawa vertex according to: 



(TTaiqi) N\Hl2\7Tb{<l2) N) 

{iTa{qi)^b{q2)N\Hl2\N) 
{N\Hl2W{qi)TTk{q2)N) 

Here a and h are the pion isospin quantum numbers and 



(3.38) 



Vab 



o j^o ^abcTc 

^Wg-lWgi 



(3.39) 



The OPEP between the pair of the first two and the third nucleon factorizes out in Eq. (|3.37p if the constants 
a," fulfill the relations 



a3 



-as , a4 



2a5 



(3.40) 
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FIG. 7; Class-Ill contributions to the 4NF. For notation, see Figs. [T] and [2] 

With this choice, the expressions for V^j^ and V^j^ from box diagrams take a particularly simple form and are 
identical to the ones arising from the cross-box contributions V^j^ and Vj)^,. We further emphasize that although 
these constraints do not completely fix the corresponding unitary transformations, the results for all class-II 
3NFs and 4NFs turn out to be independent on a^. We found that the disconnected diagrams in Fig. [5] do not 
yield 4NFs irrespective of the valies of a^. Choosing a,; as specified in Eq. (|3.40p . the last two diagrams in the 
first row in Fig. [5] lead to vanishing 4NFs while the contribution from the first two diagrams in the second row 
is given by 



Vciass-II 



2gi ai ■ qi ^4- 94 

{2F^r [9? + M2] [q,l + A/2] [q^^ + M^] 



(ti • T4 T2 • T3 - Ti • T3 T2 • T4) (fi2 • 9*4 



Tl X T2 • T4 912 X 94 • (T3 



all perm. . (3.41) 



The last two diagrams in the second row of Fig. [5] do not involve reducible topologies. The corresponding 
4NFs can, therefore, be evaluated using the Feynman graph techique with no need to consider all possible 
"time-ordered" diagrams. Using the Feynman rule 

^ [t'^,5''V • (4agi -f (4a - l)(g2 + 91)) + 2 cycl. perm.] , (3.42) 

for the ttttttNN vertex and 

[S'^'S'" ( [(91 + 92)' - + 2a [AM^ - qf - ql - ql - qf^ ) + 2 cycl. perm.] , (3.43) 

for the TTTTTTTr-vertex, where the superscripts refer to the pion isospin quantum numbers and qi denote the 
corresponding pion outgoing momenta, we obtain the following result: 



^ 9a ^2 ■ 92 o'a • 93 ^4 ■ 94 

Kciass-II + ^^2] [^-.^2 + ^2] [^^2 + ^^2] ^1 ' ^3 • X4 

■,4 



^1 ■ (93 + 9*4) - 4a(?i • (92 + 93 + 94) 



Qa • 9i 0-2 • 92 0-3 ■ 93 0-4 • 94 

-Tl • T2X3 • T4 



(91 + 92 )' + 



2(2F.)6 + M2] [g-^ + M2] [q,^ + M^] [q^ + M^] 
- 2a{AMl + q^ + qi + qi + qi) + aU perm. . (3.44) 

It is easy to see that terms proportional to a cancel out each other. Thus, while the individual contributions from 
the last two diagrams in Fig. [5] do depend on the arbitrary constant a, the a-dependence drops out completely 
in the total result. Notice further that the individual contributions of these two diagrams are given in Ref. nj 
for the choice a = 0. 

• Class-Ill contributions proportional to g\. 

This class of contributions arises from diagrams with two insertions of the Weinberg- Tomozawa vertex or 
one insertion of the ttttNNNN vertex H^2 defined in the last line of Eq. p.20p . see Fig.[71 The corresponding 
operators read: 



1 \^ 2 2 "^1 1 2 2 1 1 2 2 

2^21— H 211] H 22-^ H 22 + 21-^ H 21 rjH 22^ H 22 - -^21^-^21 ^^22 ^^22 

Xl Xl X2 1 Xl Xl Xl Xl X3 X2 



2 



trl tr2 ^ ttI tt2 ttI tt2 tt2 ttI rrl Tj2 ttI _j_ tt 

■"211=r-"22T;~-"21T^^22 77^21 -"22 T^^22 T^^21 -"21 T^^22 T^'"21 T^'"22 

^7r ^TT i^-n ^ rj-n J^tt ^tt i^-n i^ir i^ix 



FIG. 8: Class-IV contributions to the 4NF. For notation, see Figs. [T] and [2l 
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^ ffl f72 tt2 Tjl 

^ jn 2 2^ . .iJ oo , -iJ oo . o ■ 



A^ 



'22 



21 



— Unn Ho-] -ffoi 

Cf ZZ TP ^r T-i ZL TP 

Z i/TT iijr 



21 



?7 + h.c. 



1 X2 \1 \2 

TT^ ffl ffl ff2 

0^22 TP ^21 TP ^21 77 ^22 
Z i^TT i/^r iijr 



(3.45) 



The class-Ill contributions are not affected by the 77-space unitary transformations. The disconnected diagram 
is again found to produce vanishing 4NF contribution. The 4NF contributions from the last two diagrams 
derived using Eq. p.45p do not vanish individually but cancel each other. This result can also be obtained in a 
simpler way using the Feynman graph technique. In this framework, one only has to consider the second graph 
in Fig. [71 The last diagram does not appear since there is no ttttNNNN vertex of dimension = in the 
Lagrangian. Due to the four-momentum conservation at each vertex and the fact that the Weinberg- Tomozawa 
vertex contains a time derivative of the pion field, the contribution of the second Feynman diagram in Fig. [7] is 
suppressed by Q"^ jvr? and, therefore, shifted to higher orders. 

• Class-IV contributions proportional to g\ Cs,t- 

This class of contributions arises from diagrams involving four vertices iJ2i and one insertion of the leading-order 
contact interactions i?|oi Fig. [51 From Eqs. (|2.14p - (|2.17p . one obtains: 



1 A^ 

T^Hli— H^i rj — H^o -;2 ^21 - o -^2\ -TT ^21 V H^i—^; — H21 



Ai 



Ai 



-^Hli^Hl^rjHl-^—Hl^—Hl-^ - -111(^1^1121—1121—1121^1121 



HtoVH, 



Xl X2 \1 



-H. 



21 jp ^^21 JPJ2 -'-'21 



-E^ 



Et, 



' Ett 



El 



1 -j^ 12 1 1 

+ 2 -^21^^21 ^-^^40 '7-^^21 ^-f^21 - 

A^ A^ A^ A^ 
+ H2i-^H2i—H2i—H^Q—H2i 

J-JTr J^TT -C/TT -i-^TT 



1 A^ A^ A^ 

- 2-^40'7^21^-f^21^-f^21^-f^21 

3 A^ A^ 

-H^^qT] H2i^ 11211] H2i^H2i 

1 A^ A^ 

g -^21 ^-f^21 V HiO V -^21 ^^21 

1 Xl )i2 \2 \1 

Ffl ffl Wl f/l 

" o-"21T^-"21T^-"40T;~-"21Tr-"21 

Z Hjt, Hj-jt tj-rr tj-rr 



T] + h.c. (3.46) 
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(a) 



(c) 



(g) (h) (i) (j) (k) (1) 

FIG. 9: Class-IV three-nucleon diagrams which do not involve self-energy insertions. For notation, see Figs. [T] and [2l 



In addition to terms listed above one has to include contributions arising from the Ty-space unitary transformation 
defined in Eqs. ((3?24l) . ((3?25)) and acting on Hi^,^ 



V 



1 1 A^ A^ 

2 ^21 '^^21 V H2i—^H2i rj - - H^i-^H^i rj H^i—H^i rj 

1 A^ A^ A^ 1 A^ A^ A^ 

- 4 H'2i^H'2i -W-^l^Ei^li ^ ^'o + 4 Hi,— Hi,— Hi,— Hi, ^ H 



77 + h.c. (3.47) 



Here, we already adopted the values for ai^2 from Eq. p.3ip . Finally, one needs to take into account contributions 
which arise from yet unconsidered unitary transformation in Eq. (|3.24|) with S = agSg and the generator S(, 
given by 



H2i^H2, rj H^„ + H^g f] H2, H2, 



(3. 



The corresponding contributions to the effective Hamilton operator read: 



A-^ A^ A-^ A-^ 

H21 ~^H2i V H2,^H2, rj H^Q + H^Q 77 -^^21 ^-^^21 V -^21 '^^21 



)l yl yl yl 

Hli -^^l, -q HIq 7] Hi,— Hi, - ^^2\ V ^40 V Hi,— Hi, 



Notice that the second possible generator 

77 



rrl ^ ff2 ^ ttI , ttI ^ tt2 ttI 
^21 ^2 HiO ^ H2, + H2, H,o ^2 ^21 



r; = 



(3.49) 



(3.50) 



does not lead to a non-trivial transformation. The constant ag can be determined unambigously by the require- 
ment that the OPEP factorizes out in the 3NF diagrams shown in Fig. [9l The energy denominators for the 
individual graphs in this figure read: 



These terms are not shown exphcitly in Eq. 113.271 1. 
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FIG. 10: Class- V contributions to the 4NF. For notation, see Figs. [T] and [2l 

2 





— "^aW 


— ^ZN — 


~^3N 




- T/j - 


Aae 






— — 


2 


-2(1- 




— ^In = 


0. 





UJiU!'. 



3 ' 



(3.51) 

Here, i labels the pion exchanged between the first two nuclcons and the third one while j denotes the pion 
which does not interact with the third nucleon. The OPEP factorizes out in the above terms if one sets 



(3.52) 



We have verified that this choice also leads to the desired 1/wi-dependence for the class- IV 3NF contributions 
involving self-energy insertions which are not shown in Fig. [9l With the ry-space unitary transformations being 
fixed as described above only the last two diagrams in Fig. [5] yield non- vanishing 4NF contributions: 



Vciass- 



IV 



4Ct 



9\ 



ax ■ gi g3 X 0-4 ■ (712 

(2F.)4 + +M2] 



all perm.. (3.53) 



Class- V contributions proportional to g\ Gs.t- 

The class-V contributions arise from diagrams constructed from two vertices i?2ii ^'^^ Weinberg- Tomozawa 
vertex and the leading contact interaction i?|o' Fig. [TUl From Eqs. (|2.14p - (|2.17p . one obtains: 



1 A-^ 1 A-^ A^ 1 A^ A-^ 

-H^qT] H2i—H2i-^H22 + -^H^^i] H2i — H22^H2i + -j^H^i^ri H 22— H 21^ H 21 



1 A^ A^ 1 A^ A^ 1 A^ A^ 

-HIqJ] H2i^H2i—H22 + -^HIqI] H2i^H22—H2i + -^H^^V H22^H2i — H2i 



\1 X2 \2 

Wl Wl Ff2 

^ 2 1 i=r 2 1 i=r 4 i=r 2 2 

J^TT i^TT ^TT 



xi xi xi xi Xi X2 

Wl TJ"^ fj2 F/l _ F/l TT^ P/l W2 

-"211^ -"22 "FT -"40T^ -"21 -"21T^-"40T^-"21T^^22 

J^TT J^TT -C/TT J^TT J^TT i^TT 



(3.54) 
f] + h.c. . 



In addition, one has to take into account the contributions arising from the 77-space unitary transformations 
proportional to as, 4, 5 and acting on HIq-. 



"5 '7 



A^ A^ A-^ A-^ 

HiO V i?21 ^^22 ^ ^21 ~ ^40 V H21 -^^22 H21 



2 + 2^5 ) ?7 



y2 y^l ^1 y2 

H40 V H22^H2i — H21 - i?4o r] H21—H21 i?22 



FIG. 11: Class- VI contribution to the 4NF. For notation, see Figs. [T] and [D 
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a,5 ?7 



1 u-i 



El' 



Ett 



'22 



T] + h.c. 



(3.55) 



Here, we have adopted the values for a^^i from Eq. (|3.40p . Similar to the class-II forces, the resulting 4NF turns 
out to be as-independent. From the diagrams shown in Fig.[TOl only the last two in the first row generate the 
non-vanishing contribution which reads: 



Vc 



lass— V 



2Ct 



9l 



(Ji ■ gi g3 X 0-4 ■ qi2 



Ti yi T2 ■ + all perm.. 



(3.56) 



Class- VI contributions proportional to Cs,t- 

The only way the class- VI contributions can be generated is from a single disconnected diagram shown in Fig.fTT] 
The corresponding terms in the effective potential read: 



1 \2 \2 

U'^ f/2 W2 

0-"22 77. -"40 TP 22 



1^2 A2 ^2 

2 



22 j^2 ^22 V ^40 



77 -|- h.c. . 



(3.57) 



It is easy to verify that the diagram in Fig. [TT] leads to a vanishing 4NF. 
Class- VII contributions proportional to g\C^T- 

The class- VII contributions arise from diagrams involving two vertices i?2i ^^.'^ two insertions of the leading- 
order contact interactions i?4o, see Fig. [121 From Eqs. (|2.14p - (|2.17p . one obtains: 

Trf41 ^ 1 A^ 2 A^ 2 A^ 1 1 1 A^ 9 A"'^ 2 2 1 A^ 9 A^ -[^ tt2 

' = 1-] - -i?21 -f^40 H21 + -i?21 H^Q— 112111 H^Q + -H21—HIQ iJji 7?ii40 

Z LO LO UJ L LO CO" Z UJ UJ 

1 \1 

77 -f h.c. (3.58) 



Further terms arise from the T^-space unitary transformation proportional to ag, see Eq. p.48p . and acting on 
2 

40- 



HIq. Using the value for ag from Eq. p.52p . these contributions take the form 



1 



n A^ 



1 A^ 



(3.59) 



From all diagrams shown in Fig. I12| only the third, fourth and fifth ones in the second row lead to the nonvan- 
ishing 4NF: 



Vc 



Class-VII 



2(j2 9 A 0-1 X 0-2 • qi2 0-3 X 0-4 • gi2 



T2 ■ '''3 + all perm. 



(3.60) 



Class- VIII contributions from disconnected graphs with insertions of a = 2-interaction. 

Finally, we need to consider contributions involving insertions of the — 2-interactions i?|o' -^02 j ^21 ^^'^ 
iJ^o and arising from disconnected diagrams shown in Fig. 1131 The corresponding terms in the effective Hamil- 
ton operator are listed in Appendix [Aj Evaluating matrix elements of these terms we have verified that all 
disconnected diagrams in Fig. [T3| lead to vanishing contributions to the 4NF. 
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FIG. 12: Class- VII contributions to the 4NF. For notation, see Figs. [T]and[2l 



D. Discussion 



In the previous section we have worked out the leading 4NF which is given in Eqs. (|3.33p . (|3.4ip . (|3.44p . (|3.53p . (|3.56p 
and (|3.60p . Some of these contributions have, in fact, already been considered in the past. In particular, McManus 
and Riska [29] discussed a long time ago the 4NF generated by the last two diagrams in Fig. \5\ These terms were 
also studied by Robilotta |30| who, in addition, considered effects due to intermediate delta excitations and exchange 
of rho and axial-vector mesons. In that work, the ttttttNN and TTTTTTTr vertices were parametrized in terms of the 
constant ^ which plays a role similar to a in Eq. (j3.19p . Our result for the penultimate diagram in Fig. [5] shown in 
the first line of Eq. (|3.44p agrees with the one of Ref. [29] if one chooses a = 1/4 and with the one of Ref. [30] if 
one chooses 4a = ^. The result for the last diagram in Fig. [5] which is given in the last two lines of Eq. p.44p does, 
however, only agree with the one of Ref. [s^l if one sets a = ^ = 0, see Eq. (43) of that work. This disagreement can 
be traced back to the (slightly) different Feynman rules used for the TTTTTTTr vertex in that work. Notice further that 
the ^-dependence does not drop out completely in the expressions for the 4NF of Ref. All other nonvanishing 

4NF contributions discussed in the previous section result from diagrams involving reducible topologies and have, to 
the best of our knowledge, not been considered before. 



In order to test the effects of the 4NFs in few-nucleon systems explicit calculations will need to be performed. To 
get a rough idea about the size of the 4NF contributions to e.g. the a-particle binding energy, one can look at the 
strength of the corresponding r-space potentials which are given in appendix IB] The strengths of the contributions 
which do not involve contact interactions is given by 



50 keV and 



35 keV , 



(3.61) 



see also [30|- These numbers set the typical scale for the expectation values of the corresponding operators. Recently, a 
more careful estimation of the 4NF effects in ^He was carried out by Rozpedzik et al. [3l| which, however, still involves 
severe approximations in order to simplify the calculations. Expectation values of various individual long-range 4NF 
contributions (i.e. the ones which do not involve contact interactions) obtained in that work agree reasonably well 
with the above estimation. For example, for the Gaussian wave function their magnitudes lie in the range ~ 70 . . . 200 
keV. Interestingly, one observes strong cancellations among various contributions. The cancellation (to a large extent) 
between the contributions corresponding to the last two diagrams in Fig. O was already mentioned in Ref. [29|. Even 
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FIG. 13: Class- VIII contributions to the 4NF. Solid squares denote vertices of order v = 4. For the remaining notation see 
Figs. [Hand [2] 

more surprising, one observes a nearly complete cancellation for the total long-range contribution Vciass-i + Vciass-ii-^ 
For all three wave functions employed in Ref. [3l|, its expectation value does not exceed 4 keV. Further, a similar 
cancellation occurs for the class-IV and class-V contributions. For the Gaussian wave function, the corresponding 
expectation values are [sij 

(*4ffe|^Class-Iv|*4He> + (^^^e |^^Class-V l^^^e) = (14.2 kcV) X Ct - (14.9 kcV) X Ct ■ (3.62) 

Here the value for the LEG Ct should be taken in units of GeV^^. The LEG Ct has been determined from fits to 
NN phase shifts. Its value depends on the cutoffs employed and is typically of the order \Ct\ ~ 10 GeV~^ [2]. Values 
obtained from resonance saturation using various phenomenological potentials are of a similar size js^ . One therefore 
concludes that the contribution to the ''He binding energy from 4NF terms proportional to Ct is rather small. The 
most significant contribution is produced by the class- VII 4NF which is quadratic in Ct- For the Gaussian wave 
function it is found to be [3l| 

(5'4^g|Vciass-VII 

15-4^^) = (-1.1 keV) X C| . (3.63) 

For \Ct\ ~ 10 GeV~^ the 4NF is, therefore, expected to produce ~ 100 keV attraction in ^He. We, however, emphasize 
that the expectation values for the contributions involving contact interactions might be significantly overestimated 
due to approximations made in Ref. (sil for the ''He wave functions (33l |. 



IV. SUMMARY AND CONCLUSIONS 



In this paper we have analyzed the next-to-next-to-next-to-leading contribution to the nuclear Hamiltonian in chiral 
effective field theory using the method of unitary transformation. The pertinent results can be summarized as follows: 

i) We have presented a new formulation of the chiral power counting which is particularly suitable for applications 
based on algebraic approaches such as the method of unitary transformation. It allowed us to express the 



^ Expressions for and in Refs. Illl and [sill should be taken with opposite sign. 
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recursive (formal) solution of the decoupling equation in a compact form and to drastically simplify the derivation 
of the effective Hamiltonian as compared to the formalism given in Ref. [l^ . 

ii) Using this new formulation we have worked out the N'^LO contribution to the effective nuclear Hamiltonian in 
chiral EFT which involves 2N, 3N and 4N operators. The 4N operators yield the leading contribution to the 
4NF. 

iii) We employed a large class of additional rj-space unitary transformations that affect N^LO terms in the effective 
Hamiltonian. The appearance of such transformations is a new feature at this order in the chiral expansion. 
We found that certain classes of the 3NF at N'^LO cannot be renormalized for an arbitrary choice of the 
ry-space transformations. The renormalizability requirement leads to constraints for these additional unitary 
transformations but does not fix them completely. 

iv) We have worked out the leading 4NF both in momentum and configuration spaces. Enforcing the above men- 
tioned renormalizability constraints we found no umbiguity in the 4NF due to the employed 77-space unitary 
transformations. The resulting 4NF is local and depends on the pion decay constant F^^, the axial pion-nucleon 
couplingg^ and the lowest-order NN contact interaction whose strength Ct is determined in the two-nucleon 
system g^g. 

In the future, the effects of the 4NF should be tested in the four-nucleon continuum and in the spectra of light nuclei. 
Further, it is also important to study effects due to intermediate delta excitations which are expected to be significant 
[30| . Work along these lines is in progress. 
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APPENDIX A: CONTRIBUTIONS OF THE ORDER Ai = 2 VERTICES TO THE EFFECTIVE 

HAMILTON OPERATOR AT ORDER u = i 



In this appendix we consider the operator structure of the effective Hamiltonian at order v = A which results from 
higher-order vertices. As already pointed out, these terms lead to disconnected diagrams which are shown in Fig. [T51 
As a general rule, such disconnected diagrams do not yield nonvanishing contributions to the effective Hamilton 
operator in the method of unitary transformation. For the sake of completeness, we list below the complete operator 
structure of the class- VIII contributions. The relevant vertices at order = 2 are i?|g, Hq2, -f^li ^'^d 7?|g. The 
interaction together with the corresponding pion tadpole diagrams renormalize the pion mass and wave function. 
As discussed in [21], the Hq2 vertex does not show up provided the Lagrangian is formulated in terms of renormalized 
pion fields and normal ordering is applied. We will, therefore, not consider the contributions involving H^2 in what 
follows. 

Let us begin with the contributions which involve insertions of the nucleon kinetic energy 77|q . Using the values for 
ai from Eqs. (P3T|) . (1X1^)) and (^3^ one obtains 



• terms cx g\: 



1 1 1 12 1 1 1 2 \ 1 1 

V ^ r] 2 ^21^^21 ^-f^2i^-f^2i ^^20 - H2i — H2iriH2i^H2oX i?2i 



1 



H. 



A 



Ai 



A2 



Ai 



2"''K 
1 



H21 V i?20 V H21 -p- H21 - - H21 — H21 —H2i-^ H21 rj H. 



H. 



H"! Z17I \1 Ijl 

^21~-"21^-"20 -"21 



Ai 



A 1 A 



E., 



2 ^21 -^H2i^H2i—H2i 7] H. 



20 



H. 



21 



ttI ^ tt2 a2 ttI Tjl 
21^-"20'^ 21^^21 



3 A^ A^ 



.1 



8 -^2\ ;^-^2\ V H20 V H21 -^^21 



yi \2 
^Hl,—Hl,—Hl, — HUHl, 



rj + h.c. 



(A.l) 



• terms cx g\: 



V = 



A^ A^ A^ A^ 

{l + 'ia^o)H2i— H 21-^ 1122^11 - H2i—H2i^H2q\ H22 

1 A^ A^ A^ A^ 

-{I - 2a5) H2i—H22^H2i'l] H20 - H2i—H22-^H2o\ H21 



1 



A 



A2 



{l + 2a^)Hl,^Hl,—Hl2VH, 



'20 



i(l + 2^5) Hl,^Hl2^Hl,i^Hl^ 



X2 

frl tt2 rrl ^ rr2 

^21^-"20-^ -"21^-"22 



i(l - 2^5) Hl2^^Hl,^Hl, ,?i/|o - ii?2\^ff2i ^^^lo '7^^: 



Ai 



Ai 



h.c. 



(A.2) 



terms cx Cs^t: 



F72 [72 

1 xi 

i ^1 ^ ul ^ tr2 „ t72 



£^5 



40 



xi xi 

rrl ^ 772-^ ff2 i,l rrl 



40 p2- 
-^20 'f^2l'7'f^40 



1 A^ 

- H^Q TlH21^H2i7] H20 



(A.3) 



ry + h.c. 



It should be emphasized that one can construct 77-space unitary transformations with the generators depending on 
HIq. Such transformations would produce further contributions which are not considered here. 

Next consider the operators which involve one insertion of the subleading pion-nuclcon vertex iJfi • One obtains 
• terms cx 5^ : 



V 



E^ 

• terms cx gA Cs,t' 



'1 A^ A^ 
2 -^^21 ^^21 VH2i^H2i 

ttI ^ ttI ^ tt3 ^ ttI 



Ai 

" EjT 



H21 —HliTj H21 



1 1 ^ 1 1 

2 -^21 ^-^21 V H21 



\l \2 yl 

H2i-^H2i — H2i—H2i (A. 4) 

£/7r i^TT -Cj-it 



V ^r) 



rl ^ 



1.1 



H2iTrHiQ—H2i - 7; H 2i-f^ H 21 rjH^f^ - 



E^ 



E: 



El 



-El' 



40 



■q + h.c. . 



(A.5) 



Finally, the contributions involving iJ^Q read: 



xi xi 

F/l ffl 

^21 " -fJ/in „ -fJQi 



E-n 



Ett 



(A.6) 



We have verified via explicit calculations that matrix elements of all terms listed above and corresponding to 4N 
diagrams shown in Fig. 1131 vanish. 
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APPENDIX B: CONFIGURATION-SPACE REPRESENTATION OF THE FOUR-NUCLEON FORCE 



The derived expressions for the 4NF depend only on the momentum transfer of the nucleons and, therefore, take the 
local form in the configuration space: 

in' rs' rs' ri I V^n \ n ) = 6(7^ - ri ) 5{r2' - ) 8{t^ - fg ) 6{t^ - ) V^n ■ (B.l) 

Clearly, the locality of the 4NF only holds if one uses a regulator function which depends on momentum transfers 
qi as well. In this appendix we give the configuration-space representation of the 4NF for a general form of a local 
regulator. 

Depending on the topology, the 4NF is expressed in Eqs. (ESS]), (EUl), dSilll), (ESSl), (P^S]) and ([glBO)) in terms of 
different sets of momenta, namely {02, 93, (74}, {(fi, (il2^ 94} a-nd {(fi, q2, Qs, q4:}- The corresponding configuration- 
space expressions have the form: 

= / lll^ a). 

- j j ^'m^m-m- """" ^'"-f*- «' • f''^^' 

where = fi — fj. Here we have taken into account the overall factor {2tt)^S{P' — P) with P and P' being the total 
initial and final momenta of the nucleons, which is not shown explicitly in expressions of sec. IIII CI Using Eq. (|B.2[) 
we obtain: 



{ti ■ T4T2 ■ T3 ~ Tl ■ T3T2 ■ T4) V12 • V23 V43 • V23 



^Class-I - (16^f2)3 • ^12Cr4 ■ V43 



Tl X T2 • T4 V12 • V23 V23 X V43 • i?3 + Tl X T3 • T4 V43 • V23 V12 X V23 • l?2 



Tl • T4 V23 X V12 • (72 V23 X V43 • 0^3 



C/i(a;i2) C/2(a;23) Uiix^s) + all perm. 



2glMj v4 ^ 

V'Class-II — 7-7 — f^i • V12 (74 • V43 
(167rP2)3 



(Xi ■ T4T2 ■ T3 - Ti ■ T3 T2 • r4) V23 • V43 + Ti X r2 • T4 V23 X V43 • CT3 

X Ui{xi2) Ui{x23) Ui{x43) 



^ ,2^3 Tl • 1-2X3 • X4 (72 • V21 CT3 • V31 (74 • V41 CTl ■ (V31 + V41) Ui{x2l) Ui{x3i) Ui{xii) 



(16^P2) 



■ Ti • T2 X3 • X4 (1 - V? - - 2Vi • V2) CTl • Vi CT2 ■ V2 CT3 ' V3 ^4 ' V4 



87r(167rP2)3 

X / (fxoUi{xio)Ui{x2o)Ui{x3o)Ui{xiQ) + all perm. , 



Vciass-IV = /^g ^2\9 '^1 ■ (?3 X (74 • V23 

(iDTri*^)"^ 



Tl ■ T3 V12 X V23 ■ CT2 - Ti X r2 • T3 V12 • V23 

X Ui{xi2) U2{x23) W^(a;43) + aU perm. , 
Vciass-v = --^^^^ Tl X T2 • T3 (7i • V12 (73 X (74 • V23 C/i(a;i2) C/i(x23) W{xi3) + all perm 



„2 (j2 jy,j7 

Vciass-vii = — f „ "^^J T2 • 1-3 (71 X (72 • V23 (73 X CT4 • V23 W^(a;i2) J/i(x23) VF(x43) + all perm. . (B.3) 

iDTTi'^ 

Here, Xi = M-^fi, Xij = Xi ~ Xj, (V^) acts on Xi (xij) and Xij = \xij\. Further, the functions J7i_2 and W are 
defined as: 
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where J^j^2 3 denotes the corresponding regulator functions. 
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